Search for the decay </> — ► K^K 7 with the 
KLOE experiment 



The KLOE Collaboration 

F. Ambrosino c ' d , A. Antonelli a , M. Antonelli a , F. Ar chilli 11,1 , 
P. Beltrame b , G. Bencivenni a , S. Bertolucci a , C. Bini f,g , 

C. Bloise a , S. Bocchettaj' k , F. Bossi a , P. Branchini k , 
G. Capon a , T. Capussela a , F. Ceradinij' k , P. Ciambrone a , 
E. De Lucia a , A. De Santis f ' g , P. De Simone a , G. De Zorzi f ' g , 
A. Denig b , A. Di Domenico f ' g , C. Di Donate d , B. Di Micco j ' k , 
M. Dreucci a , G. Felici a , S. Fiore f ' g '*, P. Franzini f < g , C. Gatti a , 

P. Gauzzi f ' g , S. Giovannella a , E. Graziani k , G. Lanfranchi a , 
J. Lee-Franzini a ^, M. Martini a ' e , P. Massarotti c ' d , S. Meola c,d , 

S. Miscetti a , M. Moulson a , S. Miiller b , F. Murtas a , 
M. Napolitano c ' d , F. Nguyen j ' k , M. Palutan a , E. Pasqualucci g , 
A. Passeri k , V. Patera a ' e , P. Santangelo a , B. Sciascia a , 
T. Spadaro a , M. Testa f ' g , L. Tortora k , P. Valente g , 
G. Venanzoni a , R.Versaci a ' e , G. Xu a ' m 

^Laboratori Nazionali di Frascati delVINFN, Frascati, Italy. 
h Institut fur Kernphysik, Johannes Gutenberg - Universitdt Mainz, Germany. 
c Dipartimento di Scienze Fisiche dell'Universita "Federico II", Napoli, Italy 

d INFN Sezione di Napoli, Napoli, Italy 
Dipartimento di Energetica dell'Universita Sapienza di Roma, Roma, Italy. 
f Dipartimento di Fisica dell'Universita Sapienza di Roma, Roma, Italy. 
g INFN Sezione di Roma, Roma, Italy. 
h Dipartimento di Fisica dell'Universita "Tor Vergata", Roma, Italy. 
1 INFN Sezione Roma Tor Vergata, Roma, Italy. 
3 Dipartimento di Fisica dell'Universita "Roma Tre", Roma, Italy. 
k INFN Sezione Roma Tre, Roma, Italy. 
e Physics Department, State University of New York at Stony Brook, USA. 
m Institute of High Energy Physics of Academica Sinica, Beijing, China. 



Preprint submitted to Elsevier 



18 June 2009 



Abstract 



We have searched for the decay 4> ~~ ^ K K ^ 5 by detecting Kg pairs plus a photon 
and with the Ag-mesons decaying to 7r + 7r - , in a sample of about 1.5xl0 9 ^-decays 
collected by the KLOE experiment at DA$NE. The reaction proceeds through the 
intermediate states /o(980)7, ao(980)7. We find five events with 3.2 events expected 
from background processes. We obtain the upper limit: BR(<j) — ► K°K 7) < 1.9 x 
1(T 8 at 90% C.L. . 
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1 Introduction 

We have searched for the decay 4> — > K°K 7 with the KLOE experiment at 
DA$NE. This decay has never been observed. In the decay the K°K pair is 
produced with quantum numbers J PC = ++ , so that the two-kaon state is 
symmetric, and can be written in terms of K$, Kl as 



The K°K state can be in both singlet and triplet isospin state, so that the de- 
cay proceeds mainly via intermediate /o(980) (1=0) and ao(980) (1=1) scalar 
mesons, 



Thus, a measurement of the (p — > K°K 7 branching ratio (BR) can give 
information on the structure of light scalar mesons. The evaluation of the 
branching ratio depends on the way the scalar field dynamics is introduced 
and on the size of the couplings of the scalar mesons to the kaons. Moreover, 
interference between the fo and ao amplitudes may also contribute to the 
— > K°K 7 decay amplitude. Predictions found in literature are mainly based 
on the decay scheme (|2J), and give values ranging over a large interval, from 
few xlCT 9 to 1CT 7 [!T2To1l4T5][6Xrf8l9lT0] . In this respect, a similar decay into 
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charged kaons, <fi — > K + K 7, gives less information since the scalar dynamics 
is masked by the final state radiation. 

The signature of the radiative decay of the <fi meson into the state ([1]) is the 
presence of either two K$ or two Kl (and a low energy photon). KK invariant 
mass is limited by the 4> mass (1020 MeV) and by twice the kaon mass (995 
MeV). The photon energy is a function of the invariant di-kaon mass Mkk 
and is given by 

~ 2m, ■ (3) 



We search for final states with a K$Ks pair, with both Kg decaying to 7r + 7r . 
This corresponds to a reduction of the observable rate of 

l - x [BR(K S - 7T+7T-)] 2 ~ 2/9, (4) 



where the factor 1/2 accounts for the choice of the K$Ks final state among 
the two possibilities for the K°TC° state. We search for two decay vertices very 
close to the e + e~ interaction region, both having two tracks with opposite 
charge and invariant mass of the charged secondaries equal to the kaon mass. 
The invariant mass of the kaon pair must be smaller than the <fi mass. A 
low energy photon must be detected with a momentum compatible with the 
kinematics of the event. 



The expected photon energy spectrum (E-y) is shown in fig. [H corresponding 
to a KsKs invariant mass (Mkk) spectrum generated assuming only phase- 
space and radiative decay dynamics: 
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(5) 



The highest rate is expected for events with photon energy close to the maxi- 
mum value (see eq. [3]). The selection criteria have been optimized in order to 
maximize the sensitivity to the searched signal [TT] . 



2 The KLOE experiment 



DA$NE is an e + e~ collider running at a centre-of-mass energy y^i = m^, = 
1.02 GeV. The beams collide with a crossing angle of (it - 0.025) rad which 
gives a <fi meson momentum of ~15 MeV (12 MeV for a small part of the 
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Fig. 1. Expected <f> — > K°~K 'y photon energy spectrum in the 4> centre-of-mass, 
according to phase-space and radiative decay dynamics. 

dataset) towards the collider centre. The KLOE detector consists of a large- 
volume cylindrical drift chamber [12] (3.3 m length and 2 m radius), operated 
with a 90% helium - 10% isobutane gas mixture, surrounded by a sampling 
calorimeter [13] made of lead and scintillating fibers providing a solid angle 
coverage of 98%. The tracking chamber and the calorimeter are surrounded 
by a superconducting coil that provides an (approximately) axial magnetic 
field B = 0.52 T, the axis being the bisector of the external angle of the e + 
and e~ beams. The drift chamber has a momentum resolution <j{pj_)/p± ~ 
0.4%. Photon energies and arrival times are measured in the calorimeter with 
resolutions a E /E = 5.7%/v / ^(GeV) and a t = 57 ps/v^GeV) © 100 ps. The 
trigger [H] is based on the detection of at least two energy deposits in the 
calorimeter above a threshold that ranges between 50 and 150 MeV. 

For this search the full data sample acquired at the 0-peak energy by the 
KLOE experiment between years 2001 and 2005 has been used. After quality 
selection, a data sample corresponding to an integrated luminosity of 2.18 fb _1 
is available for the present analysis. The luminosity is measured with an error 
of less than 1% by counting large-angle Bhabha scattering events |15j . 

The KLOE Monte Carlo simulation (MC) package, GEANFI [16], has been 
used to provide high-statistics samples. The total number of events in each 
sample is generated with an integrated luminosity scale factor (LF) defined 
as LF = integrated luminosity (MC generation) /integrated luminosity (data). 
This scale factor is different for different <j) decay modes and different data 
taking periods. All decays have been simulated with LF=1 or 2, CP- violating 
rare decays, such as K$Kl — > 7t + 7t~7t + tt - , have been simulated with LF=10 
and 20. DA$NE operating conditions during data taking, machine parameters 
and background are adjusted in the MC on a run-by-run basis. 

Simulations of all <ft decays and CP-violating rare decays have been used to 
evaluate sources of background mimicking <fi — > K°K 7 decays. The initial- 
state radiation (ISR) process e + e~ — > ^K^K is taken into account in all the 
simulations. A further possible background source is the continuum production 
of charged pions through the process e + e~ — > 7t + 7t~7t + 7t - (7). A pure phase- 
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Fig. 2. Distributions of the two-kaon invariant mass (left) and the variable 

(right). Top: <p — > K°K°j MC. Bottom: data (dots) and background MC (grey). 
Selection cuts are shown. Background MC samples are normalized to data taking 
into account the different LF values. 

space generation with LF=10 has been used to simulate this kind of events. 
The measured cross section of this process in the <fi region [I7J has been used 
to normalize MC to data. 

To simulate the signal we generated a MC sample of 10 4 <fi — > KsKs^f — > 
7r + 7r~7r + 7r~7 events, developing a modified version of the PHOKHARA5 gen- 
erator [18]. Phokhara is a Monte Carlo event generator which simulates the 
electron-positron annihilation into hadrons plus an energetic photon from ini- 
tial state radiation (ISR) process, at the next-to-leading order (NLO) accu- 
racy. This simulation relies on general assumptions for the branching ratio 
dependence on the photon energy and on phase space, according to eq. [5j 



3 Data analysis 



We search for KsKs^y events by requiring two pairs of tracks with opposite 
charge, each pair originating from a vertex contained inside a cylindrical vol- 
ume centered on the interaction point and the beam line, 16 cm long in z and 
of 3 cm radius. We reconstruct the invariant masses m\ and m 2 , for each pair 
of tracks, using the pion mass hypothesis. The distribution of signal events in 
the plane mi, rri2 is well contained inside a circle of a few MeV radius centered 
on the Ks mass. We require the events to satisfy a 4 MeV cut on this radius. 

Once the two Ks mesons are identified, we build the invariant mass of the 
Ks pair, Mkk- In order to reject background from (p —> KsKl — > 7r + 7r~7r + 7r~ 
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Fig. 3. f2 angle between photon direction and missing momentum direction. Top: 
4> — ► K°K 7 MC. Bottom: data (dots) and background MC (grey). Selection cuts 
are shown. Background MC samples are normalized to data taking into account the 
different LF values. 

events, we retain only those events with K$Ks invariant mass below 1010 
MeV. Moreover, in order to reject badly reconstructed K$ semileptonic decays, 
which could mimic a two-pion decay, we define the missing mass m 7 using 
four-momentum conservation: 

ml = (P, - P Ksi - P Ks2 )\ (6) 



where P^, Pr S i an d Pr S2 are the 4- vectors of the <p and of the two reconstructed 
Ks- This variable must be zero in case of two Ks,l to 7t + tt~ decays, even with 
an extra photon, and is required to be zero within ±500 MeV 2 . Figure [2] shows 
the distribution of signal simulation, data and background simulation for Mkk 
and M 2 . Ninety five events survive these cuts. 

We then look for the radiative photon, by searching for a calorimeter cluster, 
without any associated track. The cluster must be in time with the observed 
event, and its position in the calorimeter must be compatible with the direction 
of the missing momentum. The angle Q between the missing momentum and 
the photon direction, derived from the cluster position, is required to be below 
0.6 rad. Figure [3] shows the Q distributions for data and MC simulations. 

We verified the correct simulation of low-energy clusters in the MC using the 
cj) — > 7t + 7t~7t° control sample: a time resolution of 520 ps is achieved for low 
energy clusters in both data and MC simulation, well reproduced by MC signal 
simulation. 

After the selection, 33 events are present in the non-rescaled MC samples 
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Fig. 4. Signal detection efficiency as a function of Mkk in the allowed range. 



is present in the MC sample for all decays. The overall efficiency for the 
signal is (24.8±0.5)%, evaluated from signal simulation. The efficiency is shown 
in fig. H] as a function of Mkk -The evaluation of this efficiency can vary 
depending on the scalar field dynamics and the size of couplings used for 
the signal simulation. To study this effect the kaon-loop approach [6] has 
been used, varying the couplings gsKK, with S=f ,a in ranges larger than 
the present spread of experimental measurements and theoretical evaluations. 
Signal efficiency has been recalculated throughout this range and compared 
with the one used in the analysis: the variation is within the statistical error. 

To normalize MC distributions to data we have estimated effective scale fac- 
tors, one for each MC sample, fitting the MC distributions for Mkk and to 
the corresponding data distributions. The results of the fits are shown in fig. [5j 
For this purpose we used a control sample made of events surviving all the 
selection requirements, except the M K k and rr£ cuts. For the simulation of all 
4> decays, the scale factor obtained from the fit is equal to 0.47 ± 0.04, while 
the MC scale factor for CP- violating rare decays and e + e~ — > 7r + 7r~7r + 7r~(7) 
decays is 0.096 ± 0.015. These factors differ from the LF scale factors by at 
most 15% and are in agreement within errors. Rescaled background events are 
listed in tab. [TJ 

The expected number of background events, once scale factors obtained from 
fits are applied, is 3.2±0.7 ; the error is the combination of the statistical error 
and the systematic error due to the fit. The number of residual events in the 
data sample is five. These numbers will be used for the upper limit evaluation. 



Background events 


Rare kaon decays 


e + e — > 7T + 7T 7T + 7T (7) 


non-rescaled 


26.0±5.1 


7.0±2.6 


rescaled by LF 


2.13±0.42 


0.70±0.26 


rescaled by fit result 


2.50±0.62 


0.67±0.27 



Table 1 

Expected background events, before and after rescaling. Errors are the combina- 
tion of both statistical and systematic contributions. Residual events obtained by 
rescaling by fit or LF are in agreement within the errors. 



for CP-violating rare decays and e + e 



7r + 7r 7r + 7r (7) decays. No event 
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Fig. 5. Data (points) and MC background fit result (solid) for Mkk (left) and m 2 
(right) distributions. 

4 Result 



The upper limit on BR(0 — > K°K 7) is evaluated as follows: 

JLdt x a(e+e- -> <p) x 1/2 x (BR{K S -> tt+tt")) 2 x e' 1 J 



/ Ldt is the integrated luminosity, 2.18 fb _1 . cr(e + e~ — > 0) is the cross section 
for production from e + e~ annihilation, equal to 3.09 /xb; this is evaluated 
from the fit to the visible cross section for the process <fi — > 777 with the data 
taken around the <j) peak, with a fractional error of 2.4%. BR(Ks — > 7r + 7r~) is 
taken as 0.692 [19]. The signal efficiency e is (24.8 ± 0.5)%. The upper limit 
on the number of events at 90% CL, S90CL, is calculated using the so-called 
"unified approach" [IT] , using the number of residual data events, and the 
residual background events reduced by la to take into account the systematic 
error on the background evaluation in a conservative way; we get S90CL = 7.49. 
Our result for the upper limit is: 

BR((f) -> K°Z°7) < 1.9 x 10~ 8 . (8) 



This is the only existing experimental result for this branching ratio. In fig. [B] 
the upper limit is compared with theoretical predictions [TP1I3141 5 6 7 8 9]fT0] . 
Ref. p] is based on the scalar kaon-loop model: in fig. [6] 2q and 4q indicate 
the prediction for a scalar qq structure or a tetraquark structure respectively; 
Ref. [2] describes this decay with a low-energy effective lagrangian. Ref. [3] 
only takes into account the fo contribution, without considering the do, which 
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Fig. 6. Comparison between the region excluded by the KLOE 90% CL upper limit 
on BR(cf) — > K°K°^f) (hatched area) and theoretical estimates. The numbers in 
abscissa are the reference numbers. 

explains the large BR value; the authors criticize the result of ref . [7] . Ref . [1] 
makes use of a chiral perturbation lagrangian which does not include explic- 
itly the scalar mesons in the intermediate state. Ref. [5115] calculates this BR 
through unitarized chiral perturbation theory based on resonant and non- 
resonant — > (/o + ao) coupling. In ref. [B] the couplings of the kaon-loop 
model are evaluated by using the SND experiment data [21]. Ref. [9] intro- 
duces a well-defined U(3) x 17(3) chiral model which incorporates ab initio 
the pseudoscalar and scalar mesons nonets. Ref. [10] repeats the calculation 
of ref. [5J using the first published KLOE data [2TJ22] . with and without the 
cr(600) contribution, and gives a large set of values. 

Following the kaon-loop model [6] we have evaluated predictions for BR(<f) — > 



K K 7) based on the gfoKK, gfo-Kw and g a oKK couplings obtained by KLOE [2T1l22f 
by fitting the data of </> decays into 7r + 7r~7, 7r 7r°7 and 77777. Allowing for varia- 
tions of the couplings within the corresponding errors, the predicted BRs vary 
between 0.4 xlO -8 and 6.8 xlO -8 , the lower values being still compatible with 
the present upper limit. 
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